arXiv:hep-ph/0211413vl 26 Nov 2002 


BIHEP-TH-2002-50 


Final State Interactions in D PP decays 

Medina Ablikim*, Dong-Sheng Du^, Mao-Zhi Yang^ 

Institute of High Energy Physics, Chinese Academy of Sciences, P.O.Box 918(4), 

Beijing 100039, China 


Abstract 

The two-body nonleptonic charmed meson decays into two pseudoscalar mesons are 
studied using one-particle-exchange method. The effects of the final state interactions 
are analyzed through the strong phases extracted from the experimental data. 

1 Introduction 

The study of the two-body nonleptonic weak decays of particles containing heavy (c, h) 
quarks appears to offer a unique opportunity to determine the basic parameters of quark 
mixing, and to investigate the mechanism of CP violation. However, the quarks are not free, 
they are bound in hadrons by strong interactions which are described by nonperturbative 
QCD. Solving the problem of nonperturbative QCD needs efforts in both experiment and 
theory. In the near future BESIII and CLEO-c detector will provide high precision data 
in charm physics including data on D meson decays, which will provide the possibility for 
understanding the physics in charm sector. 

It is interesting to study the weak decays of charmed mesons beyond the factorization 
approach [1]. In general, if a process happens in an energy scale where there are many 
resonance states, this process must be seriously affected by these resonances [2]. This is 
a highly nonperturbtive effect. Near the scale of D meson mass many resonance states 
exist. D meson decays must be affected seriously by these resonances . After weak decays 
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the final state particles rescatter into other particle states through nonperturbative strong 
interactions [2, 3]. Every D decay channel can contribute to each other through final state 
interaction (ESI). One can model this rescattering effect as one-particle exchange process 
[4, 5], namely the final state particles be scattered into other particle states by exchanging 
one resonance state existing near the mass scale of D meson. There are also other ways to 
treat the nonperturbative and ESI effects in nonleptonic D decays. One approach is that 
in which the FSIs are expressed by the phase shifts of the decaying amplitudes [6]. The 
other method is a flavor topology approach [7, 8], where the relative phases between various 
quark-diagram [9] amplitudes arise from the final state rescattering. 

The final state rescattering effects for charmed meson decays into two pions are studied 
using the one-particle-exchange method [10], where the magnitudes of hadronic couplings are 
extracted from experimental data on the measured branching fractions of resonances decays. 
In addition, a strong phase is introduced for the hadronic coupling which is important for 
obtaining the correct branching ratios in these decays. A similar analysis is applied to 
D —> PV decays [11], where P is a pseudoscalar meson, and V is vector meson. 

In the present work, we extend the study of the final state interactions in P —tttt 
decays to D ^ PP decays. The coupling constants extracted from experimental data 
are small for s-channel contribution and large for f-channel contribution. Therefore the 
s-channel contribution is numerically negligible in P —> PP decays. We safely drop the 
s-channel contribution in our discussion. In Sec. 2, we present the calculation within the 
naive factorization approach. The main scheme of one-particle-exchange method is described 
in Sec. 3. We give the numerical calculations and discussions in Sec. 4. The final section is 
reserved for summary. 


2 Calculations in the factorization approach 


The charmed meson decay can be described by the low energy effective Hamiltonian[12] 


n 


72 


E (CiQf + C2QI) 


( 1 ) 


where Ci and C 2 are the Wilson coefficients at rric scale, Vq is the product of Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements and defined as 


Vq = VuqV*q, 


( 2 ) 
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and the current-current operators are given by 


Ql = {uq)v-A{qc)v-A, Q\ = {Ua,q(3)v-A{qi3Ca)v-A- (3) 

We do not consider the contributions of QCD and electroweak penguin operators in the 
decays of D —> PP because their contributions are negligible in D decays. QCD factorization 
approach [13] is inapplicable to these decay modes, as the charmed meson is not heavy 
enough. The values of Ci and C 2 at rric scale are taken to be [12] 

Cl = 1.216, C 2 = -0.415. 

In the naive factorization approach, the decay amplitude can be generally factorized into 
a product of two current matrix elements and can be obtained from (1) 

A(D* ir+ir“) = V,^ (o, + 02) - ml) F°''(ml), 

A(D° ir+ir-) = ^ V„j V/j o, iA(ml - ml) F^^iml), 

A(D° wV) = Kj vs <12 - ml) f 

A(D* ^ V2 [oi - m^) F°“{rnl) 

+ <12 ifKlmo - ml) F°*(mK)] , 

A(D» ^ ^ V'S Oi tfAml - mi-) F°“(ml), 

A(D<‘ ^ A'“7r») = V'S aj 2/i'(m|, - mi) F’^l’ni), 

A(D° ^ KA-) = ^ 14. Vj «l iMmi, - mi) F'^’imi), 

A(D+ K+A) = v^. K; ai ifidml - ml) F^’Aml), 

A(D+ ^ A'»2r+) = ^ 14. Vj «2 iMmi, - ml) 

4(B» ^ = ^k i/. 

/1(D+ ^ A+A") = ^ V4. V; Oi tfK(ml - mi) F^'^lmi), 

A(D« A+A-) = ^ V4. V; Oi 2/K(mJ, - m?,) F^'^lmi), 

A{D^ K^k^) = 0, (4) 

where the parameters Oi and 02 are dehned as [8] 

ai = Cl -I- C 2 , a2 = C 2 + Cl -f , (5) 
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with the color number Nc = 3, and x is the phenomenological parameter which takes into 
account nonfactorizable correction. For dependence of the form factors, we take the BSW 
model [1], i.e., the monopole dominance assumption: 




m 

1 - g2/m2 • 


( 6 ) 


where m* is the relevant pole mass. 

The decay width of a D meson at rest decaying into PP is 

r(D^PP) = -!-|yl(D^PP)|2M (7) 

oTT TTL 

where |p | is the 3-momentum of each hnal meson. The corresponding branching ratio is 


BriD PP) 


V{D PP) 


T 


tot 


( 8 ) 


A comparison of the branching ratios of the naive factorization result with the experimental 


Table 1: The branching ratios of Zf ^ PP obtained in the naive factorization approach and 
compared with the experimental results. 


Br (Theory) 

X = -8.6 X 10-2 


Decay mode 

TT+TT^ 
7r’''7r- 
7r°7r° 
D+ K^n+ 
K_-n+ 
DO ro^o 

po K+tt- 
D+ K+tt° 
D+ K^Tr+ 
po 

K+k^ 

K+K- 

po ^ 


Br (Theory) 
X = 0 
3.1 X 10-^ 
2.48 X 10-3 
9.98 X 10-s 
1.20 X 10-3 
4.81 X 10-2 
2.93 X 10-® 
1.88 X 10-3 
2.4 X 10-3 
3.86 X 10-8 

7.58 X 10-3* 
9.15 X 10-3 

3.59 X 10-3 
0 


2.71 X 10-3 
2.65 X 10-3 

1.39 X 10-3 
9.98 X 10-2 
5.14 X 10-2 
4.1 X 10-3 
2.01 X 10-3 
2.56 X 10-3 

5.39 X 10-6 
1.06 X 10-6 
9.76 X 10-3 
3.83 X 10-3 

0 


Br (Experiment) 


(2.5 

± 

o 

X 

o 

1 

(1.43 

± 

0.07) X 10 

(8.4 

± 

2.2) X 10~ 

(2.77 

± 

0.18) X 10' 

(3.80 

± 

0.09) X 10 

(2.28 

± 

0.22) X 10' 

(1.48 

± 

0.21) X 10' 


(5.8 ±0.6) X 10-3 
(4.12 ±0.14) X 10-3 
(7.1 ± 1.9) X 10-3 


data is presented in Table 1. The second column gives the pure factorization result, where 
the nonfactorization effect is zero, while the third column represents the branching ratio 
with small nonfactorization correction. One can notice that the results are not in agreement 
with the experimental data. For doubly Cabibbo-suppressed decay modes, the experimental 
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measurements of their decay rates are unavailable, except for the channel . We 

shall predict their theoretical branching ratios in section 4. The ratio for Cabibbo-suppressed 
decay mode vanishes in the naive factorization approach. This decay seems to 

be induced through hnal state rescattering. 

3 The one particle exchange method for FSI 

As we have seen above, the experimental results for the branching ratios are mostly in 
disagreement with the calculation from the naive factorization approach. The reason is that 
the physical picture of naive factorization is too simple, nonperturbative strong interactions 
are restricted in single hadrons, or between the initial and hnal hadrons which share the same 
spectator quark. If the mass of the initial particle is large, such as the case of B meson decay, 
the effect of nonperturbative strong interactions between the hnal hadrons most probably is 
small because the momentum transfer is large. However, in the case of D meson, its mass is 
not so large. The energy scale of D decays is not very high. Nonperturbative ehects may give 
large contribution. Because there exist many resonances near the mass scale of D meson, it 
is possible that nonperturbative interactions propagate through these resonance states, such 
as, A:*(892), A:*(1430), /o( 1710), p(770), 0(1020) etc. 

The diagrams of these nonperturbative rescattering ehects can be depicted in Figs.l and 2. 
The hrst part D P 1 P 2 or D ^ V 1 V 2 represents the direct decay where the decay amplitudes 
can be obtained by using naive factorization method. The second part represents rescattering 
process where the ehective hadronic couplings are needed in numerical calculation, which 
can be extracted from experimental data on the relevant resonance decays. 

P3 


Pi 

Figure 1: s-channel contributions to final-state interactions in 11 ^ PP due to one particle 
exchange. 

Fig.l is the s-channel contributions to the hnal state interactions. Here Pi and P 2 are the 
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intermediate pseudoscalar mesons. The resonance state has the quantum number = O’*"’*' 
derived from the hnal state particles P 3 and P 4 . From Particle Data Group [14], one can only 
choose /o(1710) as the resonance state which evaluates the s-channel contribution. However, 
the coupling of /o(1710) with two hnal mesons P 3 and P 4 is too small [10], we drop the 
s-channel contribution in the numerical calculation. 



P3 


P4 



(b) 


Figure 2: t-channel contributions to final-state interactions in P ^ PP due to one particle 
exchange, (a) Exchange a single vector meson, (b) Exchange a single pseudoscalar meson 

Fig.2 shows the f-channel contribution to the hnal state interactions. Pi, P 2 and Vi, 
V 2 are the intermediate states. They rescatter into the hnal state P 3 P 4 by exchanging one 
resonance state V or P. In this paper the intermediate states are treated to be on their mass 
shell, because their oh-shell contribution can be attributed to the quark level. We assume 
the on-shell contribution dominates in the hnal state interactions. The exchanged resonances 
are treated as a virtual particle. Their propagators are taken as Breit-Wigner form 

- 2 V r ’ (9) 

4 - imi tot 

where Vtot is the total decay width of the exchanged resonance. To the lowest order, the 
ehective couplings of /o to PP and VV can be taken as the form 

Pi = ( 10 ) 

Pi = f, ( 11 ) 

where 0 is the pseudoscalar held, the vector held. Then the decay amplitudes of /o —PP 
and VV are 

TfPP = 9fPP: (12) 
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Figure 3: The effective coupling vertex on the hadronic level 

Tfvv = (13) 

The coupling constants gfpp and g/vv can be extracted from the measured branching frac¬ 
tions of /o — PP and VV decays, respectively [14], Because /o — VV decays have not 
been detected yet, we assume their couplings are small. We do not consider the intermediate 
vector meson contribntions of s-channel in this paper. 

For the t-channel contribution, the concerned effective vertex is VPP, which can be 
related to the V decay amplitnde. Explicitly the amplitnde of E — PP can be written as 

Typp = gvpp e ■ (pi — P2), (14) 

where pi and p2 are the four-momentum of the two pseudoscalars, respectively. To extract 
gjpp and gvpp from experiment, one should square eqs.(12) and (14) to get the decay widths 

r(f^pp) = 

r(V^ ^ PP) = iT I P U - 2 m? - 2 m? + (15) 

3 Stt [ niy J rriy 

where mi and m 2 are the masses of the two hnal particles PP, respectively, and | p | is 
the momentnm of one of the hnal particle P in the rest frame of V or /. From the above 
equations, one can see that only the magnitudes of the effective couplings | gfpp \ and 
I gvpp I can be extracted from experiment. If there is any phase factor for the effective 
coupling, it would be dropped. Actually it is quite possible that there are imaginary phases 
for the effective couplings. As an example, let us see the effective coupling of gK*Kn shown 
in Fig.3, which is relevant to the process K* Kn. On the quark level, the effective vertex 
can be depicted in Fig.4, which should be controlled by nonperturbative QCD. From this 
hgure one can see that it is reasonable that a strong phase could appear in the effective 
coupling, which is contribnted by strong interactions. Therefore we can introdnce a strong 
phase for each hadronic effective coupling. In the succeeding part of this paper, the symbol 
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Figure 4: The effective coupling vertex on the quark level 


g will only be used to represent the magnitude of the relevant effective coupling. The total 
one should be ge^^, where 9 is the strong phase coming from Fig.4. For example, the effective 
couplings will be written in the form of gfppe^^^^^ and 

The decay amplitude of the s-channel final state interactions can be calculated from Fig.l 

i gi 92 

k2 _ _|_ imVtot ’ 

(16) 

where pi and p 2 represent the four-momenta of the pseudoscalar Pi and P 2 , the amplitude 
A{D —> P 1 P 2 ) is the direct decay amplitude. The effective coupling constants gi and g 2 
should be gjpp or gvpp which can be obtained by comparing eq.(15) with experimental 
data. By performing integrals, we obtain 


Af = - 


d^pi 


d?p2 


{2Tif2EiJ {2Tif2E2 


(27r) V(pp - Pi - P 2 ) A{D ^ P 1 P 2 ) 


A^^^ = 


1 i Q^ Qi 


Svrmp'^"' ' " “ k‘^ — m? + imTtot 

The f-channel contribution via exchanging a vector meson (Fig.2(a) ) is 


(17) 


.FSI 

^t,V 


1 

2 


d^pi 


d^P2 


{2t^Y2EiJ (27r)32P2 
xpi ca ■ (pi +P3)- 


(27r) V(pp - Pi - P2)A{D ^ P 1 P 2 ) 


i e 


i{9i+92) 


F{k) P2 eA ■ (P2+P4), (18) 


k 2 _ jyj ^2 _|_ tot 

where E{kY = (A^ — m^)/(A^ — is the form factor which is introduced to compensate 
the off-shell effect of the exchanged particle at the vertices [15]. We choose the lightest 
resonance state as the exchanged particle that gives rise to the largest contribution to the 
decay amplitude. 

We furthermore have 


aFSI _ 
^t,V — 


P d{cos9) 
'-1 IGTirrip 


\pi\A{D ^ P 1 P 2 ) gi 


i e 


i(6»i+02) 


k2 _ trrp + imVtot 


E{ky P2 P, 


(19) 
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where 


H = - m^ij-^{ml + ml + ml + ml) + {\pi\\p^\ + \p2\\p3\)cos6 + EiE4 + E2E3 




The t-channel contribution by exchanging a pseudoscalar meson (Fig. 2 (b) )is 


aFSI 

^t,p 


^ 1 ( 91 + 62 ) 


and we obtain 


4fy = r «|p- 


J-i W^rriD — w? + imV 


j p'i{ 9 i+ 92 ) 

X g, g2 E{ef {-H, + H 2 ), 


where 


n 

Hi = AimvJvA^D + ^ 2)211 - m 3 - ml) 


—^(.Fi.F3 - IpiIIpsI cos 6*)(ElE 4 + |pi||p 4 | cos6') 
mi 

-^(^2^4 - \P2\\P4\ COs6')(E2E3 + \p2\\pz\ COS 6) 


2-^(m^ - mf - m2)(EiE3 - IpiUpsI cos6')(E2E4 - IP2IIP4I cos6') , ( 23 ) 


2mim2 


H2 = ^ ^ 2^3 + IP2IIP3I COS^ - ;^(m|, - mi - m2)(EiE3 - |pi||p3| cos^) 

(m£) + m 2 ) [ 2mi 

E1E4 + |pi||p4| cos^ - ^^(m^ - m^ - m2)(E2E4 - IP2IIP4I cos6') , ( 24 ) 

and X represents the relevant direct decay amplitude of D decaying to the intermediate 
vector pair Vi and V2 divided by (Vi|(V^ — A)^|0)(V"2|(ld — A)^\D), 

A{D ^ lhf 4 ) 

- (lh|(f^-dl),| 0 )(\/ 2 |(V^-dl)A‘|E)- 

4 Numerical calculation and discussion 

In order to calculate FSI contribution of D decays to tttt, Ktt and KK one needs to analyze 
which channel can rescatter into the hnal states. The rescattering processes are D tttt — 
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Figure 5: Intermediate states in rescattering process for D ^ tt-k decays 

TTTT, D —KK Tin, D ^ pp ^ Tin, D K*K* nn for D ^ titt decays; D —> 
Kn Kn, D — K*p Kn for D —Kn channels; D nn KK, D —KK — KK, 
D ^ Tip ^ KK, D ^ pp ^ KK, D K*K* KK, D ^ pcj) ^ KK for D KK 
decays and pictorially shown in Fig. 5, Fig. 6 and Fig. 7. These rescattering processes 
give the largest contributions, because the intermediate states have the largest couplings 
with the hnal states and the masses of exchanged meson are small which give the largest 
t-channel contributions. When we calculate the contribution of each diagram in Figs. 5 ~ 
7 via eqs.(19) and (22), we should, at hrst, consider all the possible isospin structure for 
each diagram and draw all the possible sub-diagrams on the quark level. Secondly, write 
down the isospin factor for each sub-diagram. For example, the uu component in one hnal 
meson 7r° contributes an isospin factor and the dd component contributes For 

the intermediate state 7r°, the factor and —should be dropped, otherwise, the isospin 
relation between different channels would be violated [10]. Third, sum the contributions of 
all the possible sub-diagrams on the quark level to get the isospin factor for each diagram 
on the hadronic level. 

The FSI contributions of the Cabibbo suppressed decays D~^ and 7r°7r° 

depend on the couplings and phases and 5 fp,r 7 re*^^’^"^respectively, while 
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(d) 


(f) 


Figure 6: Intermediate states in rescattering process for D ^ Ktt decays 




K 


K 



K 

K 


Figure 7: Intermediate states in rescattering process for D KK decays 
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7r~^7r~ depends on both of them. In the Cabibbo favored decays K^7r~^, 

K~7r^, and the doubly Cabibbo suppressed decays K^'k~, K~^n^, 

D~^ —^ K^71~^ and the branching ratios, including both the direct decay and 

the rescattering effect, depend not only on and but also gpKK^^^^^^■ 

For the Cabibbo suppressed decay modes K^K~ and K^K^, 

the FSI effects come from gK*K-Ke^^^*^^, gpKKe'^^^^^, and gK*Kr}e^^^*^'^ ■ 

In the numerical calculation, we use the input parameters: 1) the decay constants, /^r = 
0.133 GeV, fK = 0.162 GeV, fp = 0.2 GeV, /k* = 0.221 GeV, = 0.233 GeV; 2) the form 
factors, F^^(O) = 0.692, F^^{0) = 0.762, ylf^*(0) = 0.880, 7lf^*(0) = 1.147, 7lf^(0) = 
0.775, A2^{0) = 0.923 [1]. Except for the decay constants, the values of the form factors 
have not been known exactly yet. We therefore have to take them from model-dependent 
calculations. The parameter A in the off-shellness compensating function F{k‘^) introduced in 
eq.(18) takes the value of 0.513GeG, which is in the mass ranges of the hnal state mesons. In 
order to get the branching ratios which include both the direct decays and the rescattering 
effects, we use eq.(15) and the center values of the measured decay width of K* Kn, 
p —TTTT and (f) KK [14] to obtain gK*Kn = 4.59, = 6.0 and g^^KK = 5.77. We 

take gpKK = with the ss suppression parameter A=0.28 [16]. Since there is no data 

existing for K* Kg (the mass of K* is not large enough to decay into Kg), we estimate 
the value of the strong coupling gK*Krj — 3.5 by comparing gp^^^ with gx^K-Kt and considering 
SU{3) flavor symmetry with 20% ~ 30% violation. When the nonfactorizable parameter y 
is not taken into account, we can not reproduce the experimental data for all the D PP 
decays simultaneously. So we need to keep it as a phenomenological parameter. By taking 
X = —8.6 X 10“^, the experimental data of all the detected D —> PP decays can be well 
accommodated within the experimental errors. 

The strong phases of the effective hadronic couplings 9k*k-ki 9pKK, 9 ^kk and 9x*K'q 
can not be known from direct experimental measurement or from nonperturbative calcula¬ 
tions, because there are no any such kind of computations yet. The only information is that 
the values of these phases should not differ too much, according to SU (3) flavor symmetry. 
We fit the experimental data to get the values for these phase parameters, and find that it is 
possible to reproduce the experimental data of these D decays with small SU{3) flavor sym¬ 
metry violation. To show this situation, in which the experimental data are accommodated. 
Table 2 gives the numerical results of the branching ratios at 9k*ktt = 53.9°, 9pT,T, = 57.3°, 
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Table 2: The branching ratios oi D ^ PP. 


Decay mode 

Factorization 

Factorization -|- FSI 

Experiment 


2.71 X IQ-o 

1.8 X 10-4 

(2.5 ±0.7) X 10-4 

1 

+ 

T 

o 

2.65 X IQ-o 

1.49 X 10-4 

(1.43 ±0.07) X 10-4 

7r°7r° 

1.39 X 10-® 

1.06 X 10-4 

(8.4 ±2.2) X 10-4 

D+ K^t^+ 

9.98 X 10-2 

2.95 X 10-2 

(2.77 ±0.18) X 10-2 

DO ^ K-TT+ 

5.14 X 10-2 

3.72 X 10-2 

(3.80 ±0.09) X 10-2 

DO ^ 

4.1 X 10-4 

2.09 X 10-2 

(2.28 ±0.22) X 10-2 

DO ^ K+tt- 

2.01 X 10-4 

1.41 X 10-4 

(1.48 ±0.21) X 10-4 

D+ D+TfO 

2.56 X 10-4 

2.96 X 10-4 

- 

D+ ^ A:07r+ 

5.39 X lO-o 

7.56 X 10-4 

- 

DO ^ ATOttO 

1.06 X lO-o 

2.84 X 10-4 

- 

D+ ^ K+K^ 

9.76 X lO-o 

6.4 X 10-4 

(5.8 ±0.6) X 10-4 

DO ^ K+K- 

3.83 X lO-o 

4.0 X 10-4 

(4.12 ±0.14) X 10-4 

DO ^ K^RO 

0 

5.73 X 10-4 

(7.1 ±1.9) X 10-4 


dpKK = 71.8°, 9^kk = 58.7 and 0K*Kri = 65°, with a small SU{3) symmetry breaking effects. 
Column ‘Factorization’ is for the branching ratio predicted in naive factorization approach, 
where the nonfactorizable correction is small. We find that the data of H > PP cannot 
be accommodated without including the contribution of the nonfactorizable effect. Column 
‘Factorization + FSF is for the branching ratio of naive factorization including the final 
state interaction. The contributions of final state rescattering effects are large, which can 
improve the predictions of naive factorization to be consistent with the experimental data. 
The strong phases introduced for the effective hadronic couplings gK*Kn, gpmr, gpKK, gtpKK 
and gK*Kp are important for explaining the experimental data, otherwise, it is quite difficult 
to get the correct results for these decay modes at the same time by varying other input 
parameters. 

As we have mentioned earlier, the branching ratios have not been detected in experiment 
for doubly Cabibbo-suppressed decay modes D~^ K^n~^ and K^n^, 

except K^n~ yet. In our method, they are all predicted to be at order of 0(10“^^). 

To conclude this section, we shall give some comments. There are some free parameters, 
such as the D decay form factors which have not been well determined in experiment yet. 
They need to be measured from leptonic and semileptonic decays of D mesons, which are 
quite possible in CLEO-C program in the near future. The other input parameters that may 
cause uncertainties are the shape of the off-shell compensating function F{k‘^) and the non- 
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factorizable parameter x, which are needed to be studied by some nonperturbative methods 
based on QCD in the future. Certainly to completely understand hnal state interactions, 
more experimental data and more theoretical works are needed. 

Summary We have studied two-body nonleptonic charmed meson decays into two pseu¬ 
doscalar mesons. The total decay amplitude includes both direct weak decays and final state 
rescattering effects. The direct weak decays are calculated in factorization approach, and the 
hnal state interaction effects are studied in one-particle-exchange method. The prediction 
of naive factorization is far from the experimental data. After including the contribution of 
hnal state interaction, as well as the nonfactorizable correction, the theoretical predictions 
can accommodate the experimental data within experimental errors, where the strong phases 
of the ehective couplings are quite necessary to reproduce experimental data. The branching 
ratios are predicted for the three doubly Cabibbo-suppressed decay modes. 
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